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NA44 uses a 512 channel Si pad array covering 1.5 < rj < 3.3 to study charged hadron 
production in 158 A GeV Pb+Pb collisions at the CERN SPS. We apply a multiresolution 
analysis, based on a Discrete Wavelet Transformation, to probe the texture of particle 
distributions event-by-event, by simultaneous localization of features in space and scale. 
Scanning a broad range of multiplicities, we look for a possible critical behaviour in 
the power spectra of local density fluctuations. The data are compared with detailed 
simulations of detector response, using heavy ion event generators, and with a reference 
sample created via event mixing. An upper limit is set on the probability and magnitude 
of dynamical fluctuations. 



1. Introduction 

The deconfmement /chiral symmetry restoration phase transition in ultrarelativistic 
heavy ion collisions is inherently a multiparticle phenomenon. Therefore, event-by-event 
analysis of multiparticle hadronic observables is of paramount interest. We carry out a 
texture, or local fluctuation, analysis to determine the correlation/fluctuation content of 
single events, and analyze the scale composition of the correlations. 

The idea to look at particle distributions in rapidity y to search for critical behaviour was 
proposed |l]Q] based upon a Ginzburg-Landau type of multihadron production theory [|l|, 
where a random hadronic field <fi(y) plays the role of an order parameter in a hadronization 
transition. Enhanced large scale correlations of hadrons in y at the phase transition 
would signal critical fluctuations in the order parameter. Stephanov and coworkers 
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Figure 1. Power spectra : O ~~ true events, 
A - mixed events. 



Figure 2. Confidence coefficient clS 8l func- 
tion of the fluctuation strength. 



indicated a second order QCD phase transition point which should exist under certain 
initial conditions, within the reach of today's experiments. 

In our work, a power spectrum analysis of event texture in pseudorapidity rj and az- 
imuthal angle (, based on a Discrete Wavelet Transformation (DWT)Q, is performed on 
a number of large event ensembles sampled according to their multiplicity, thereby study- 
ing the impact parameter dependence of the observables. D WT quantifies contributions 
of different ( and rj scales into the event's overall texture, thus testing the possible large 
scale enhancement. 



2. Experiment and Results 

Ionization energy loss of charged particles was measured in 512 silicon pads, with radial 
granularity of 16 and azimuthal granularity of 32. The pads were read out by AMPLEX 
[[J chips, one chip per sector, (^-electrons, produced by the Pb beam traversing the target, 
were swept away to one side by a dipole magnetic field (< 1.6 Tl). Only the 5-electron-free 
side was used in the analysis. Empty target runs were used to measure the background. 
Cross-talk in the detector was evaluated off-line. 

Track density in an individual event p((, if) is represented by a 2D array of the calibrated 
digitized amplitudes of the channels ( an amplitude array). This is expanded into a basis 
of Haar wavelet^ functions ^^^(C, rj), orthogonal with respect to scale fineness m and 
location (i, j), with the A index numbering three modes of texture direction sensitivity - 
C, r], and (rj - for azimuthal, pseudorapidity and diagonal. f\ With our binned detector, 
the act of data taking is the first stage of the Haar wavelet transformation (not true for 
any other wavelet). 

From the squared expansion coefficients a power spectrum is formed: 



(1) 



2 Technically, the 2D basis ^ itj (C,v) = 2 m * A (2 m C 
the following way: ^ = ip(0<p(rj), ^ = 0(C) V> fa), * C '' 



function tp(x) — {+1 for < x < 
for < x < 1 and otherwise. 



-1 for ± 



z,2 m ?7 — j) is constructed from ID functions in 
: ip(Qil>(r)), where for any variable x, the wavelet 
< x < 1;0 otherwise}, and the scaling function 4>{x) = 1 
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Figure 3. Multiplicity dependence of the tex- 
ture correlation. O ~ tne NA44 data, • - RQMD. 
The boxes show the systematic errors vertically 
and the boundaries of the multiplicity bins hori- 
zontally; the statistical errors are indicated by the 
vertical bars on the points. The rows correspond 
to the scale fineness m, the columns - to the di- 
rectional mode A (discussed in the text). 



Figure 4. Coarse scale T] tex- 
ture correlation in the NA44 data, 
shown by O (from the top right 
plot of Figure |3|), is compared 
with that from the multifireball 
event generator for three different 
fireball sizes. Detector response is 
simulated. 



We used WAILI |7j software to obtain the wavelet expansions. More experimental and 
algebraic details can be found in || and the references therein. 

Figure [1] shows that the power spectra in central Pb+Pb collisions are indeed enhanced 
on the coarse scale. It is necessary to eliminate "trivial" and detector related contributions 
to this enhancement. Event mixing is done by taking different channels from different 
events. In order to reproduce the electronics cross-talk effects in the mixed event sample, 
mixing is done sector-wise, i.e. the sectors constitute the subevents subjected to the 
event number scrambling. The mixed events preserve the texture associated with detector 
position offsets, the inherent dN/ dr] shape, cross-talk and dead channels. This is static 
texture as it reproduces its pattern event after event; we are interested in dynamic texture. 
We reduce sources of static texture in the power spectra by empty target subtraction and 
by subtraction of the mixed events power spectra, thus obtaining the texture correlation 
P x (m) true — P X ( m ) mix- This quantity, normalized to the RMS fluctuation of P x (m) mix , is 
used to characterize the relative strength of local fluctuations in an event. Its distribution 
for different A is plotted on Figure ^| in an integral way, i.e. as an a(x) graph where for 
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every x, a is the fraction of the distribution above x. The confidence level with which 
local fluctuations of a strength x can be excluded is then 1 — a. Fluctuations greater 
than 3 x RMS mix are excluded in the azimuthal and pseudorapidity modes with 90% and 
95% confidence, respectively. The monotonic fall of the curve is consistent with absence 
of abnormal subsamples in the data. 

Figure ^ summarizes texture correlation data for different multiplicity events. The sys- 
tematic errors were evaluated by removing the Pb target and switching magnetic field 
polarity to expose the analyzed side of the detector to ^-electrons, while minimizing nu- 
clear interactions. Correlations (i.e. deviations of P x {m) trU e from P x (m) m i X ) in such 
events measure the remaining systematic uncertainties. Thus, this component of the sys- 
tematic error is signed, and the systematic errors are asymmetric. The other component 
(significant only on the coarsest scale) is the uncertainty of our knowledge of the beam's 
geometrical cross-section. 



3. Discussion and Conclusions 

A Monte Carlo simulation was performed, including known background texture effects 
and the subtraction of mixed simulated events. The finite beam size effect, irreducible by 
the subtraction methods, is thereby taken into account. 

The finite beam size effect explains the rise of RQMD points with dN / drj in Fig. [| 
Finally, the sensitivity of the method is evaluated (see Figure in MC by generating 
multiple fireball events with given mean number of particles (charged and neutral), N p , 
per fireball. The total px of each fireball is 0; its total pz is chosen to simulate longitudinal 
flow by Lorentz-boosting the fireballs along the Z direction, keeping the total p of an event 
at 0. By varying N p per fireball, one varies "grain coarseness" of the event texture in t]. 
The data are consistent with clustering of N p per fireball below 50. 

This novel method of event-by-event analysis, applied to the SPS PbPb data, does not 
reveal any evidence of critical phenomena. The authors thank N.Antoniou, I.Dremin, 
E.Shuryak, M.Stephanov, and T.Trainor for illuminating discussions. This work was 
supported by the Austrian Fonds zur Forderung der Wissenschaftlichen Forschung; the 
Science Research Council of Denmark; the Japanese Society for the Promotion of Science; 
the Ministry of Education, Science and Culture, Japan; the Science Research Council of 
Sweden; the US Department of Energy and the National Science Foundation. 
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